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Abstract 

d 

injector diameter 

A theoretical and experimental investigation 
of particle-laden, weakly swirling, turbulent free 

d P 

particle diameter 

jets was conducted. Glass particles, having a 
Sauter mean diameter of 39 pm with a standard 

f 

mixture fraction 

deviation of 15 pm, were used In the study. A 
single loading ratio (the mass flow rate of par- 

G 

particle mass flux 

ticles per unit mass flow rate of air) of 0.2 was 
used in the experiments. Measurements are reported 

k 

turbulence kinetic energy 

for three swirl numbers, ranging from 0 to 0.3. 
The measurements included mean and fluctuating 

*-e 

dissipation length scale of eddy 

velocities of both phases, and particle mass flux 
distributions. Measurements were compared with 

m p 

particle mass 

predictions from three types of multiphase flow 
analysis, as follows: (1) locally homogeneous 

n 

number of particle groups 

flow (LHF) , where slip between the phases was 

"1 

number of particles per unit' time In 

neglected; (2) deterministic separated flow 
(DSF), where slip was considered but effects of 

group 1 

turbulence/particle interactions were neglected; 
and (3) stochastic separated flow (SSF), where 

P 

pressure 

effects of both interphase slip and turbulence/ 
particle interactions were considered using random 

R 

injection tube radius 

sampling for turbulence properties in conjunction 
with random-walk computations for particle motion. 

Re 

Reynolds number 

For the particle-laden jets, the LHF and DSF 
models did not provide very satisfactory predic- 

r 

radial distance 

tions. The LHF model generally overestimated the 
rate of decay of particle mean axial and angular 

S * 

source term 

velocities with streamwlse distance, due to the 
neglect of particle Inertia. LHF model predictions 

s P<t> 

source term due to particles 

of particle mass flux also showed poor agreement 
with measurements due to the assumption of no-slip 

t 

time 

between phases. The DSF model also performed 
quite poorly for predictions of particle mass flux. 

te 

eddy lifetime 

because turbulent dispersion of the particles was 
neglected. The SSF model, which accounts for both 

u 

axial velocity 

particle inertia and turbulent dispersion of the 
particles, yielded reasonably good predictions 

V J 

volume of computational cell j 

throughout the flow field for the particle-laden 
jets. 

V 

radial velocity 

Nomenclature 

w 

angular velocity 

a acceleration of gravity 

X 

axial distance 

Co drag coefficient 

Xp 

particle position vector 

Cj parameters in turbulence model 

dXp 

relative path length of particles in an 
eddy 
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Atp 

particle residence time In an eddy 
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e rate of dissipation of turbulence kinetic 

energy ' " . ” v . 

nt turbulent viscosity 

p density 

a, turbulent Prandtl /Schmidt number 

4> generic property 

Subscripts : 

c centerline quantity 

m maximum value 

o injector exit condition 

p particle property 

<= ambient condition 

Superscripts : 

( )' root-mean-square fluctuating quantity 

(-) time averaged value 

(*) vector quantity 

(~) Favre averaged value 


Introduction 

Multiphase flows are extremely important in 
many practical devices. In a gas turbine engine, 
for example, fuel is sprayed into a highly tur- 
bulent recirculating flowfield, where it evapor- 
ates and takes part in chemical reactions. 
Particle-laden flows are being studied as a step 
toward a better understanding of complex multi- 
phase flows. Particle-laden flows allow the study 
of interactions between the continuous and dis- 
persed phases without interference from vaporiza- 
tion or combustion effects. The objective of the 
present investigation was to extend previous work 
on free particle-laden jets to consider effects of 
swirl, which is an important aspect of many prac- 
tical sprays. The performance of typical two-phase 
flow models was of particular interest. Earlier 
theoretical and experimental investigations of 
particle-laden jets have recently been reviewed; 1 
therefore, the present discussion of past work 
will be brief. 

Earlier investigations involving analysis and 
experiments for nonswirling particle-laden jets 
are reported in Refs. 2 to 5. A recent study of 
nonswirling particle-laden jets is reported by 
Mostafa, et al. 8 Three methods of analyzing 
multiphase flows were considered in these studies, 
as follows: (1) a locally homogeneous flow- (LHF) 

model, where properties of both phases were taken 
to be the same; (2) a deterministic separated flow 
(DSF) model, where finite interphase transport was 
considered but effects of turbulence on particle 
motion were ignored; and (3) a stochastic separ- 
ated flow (SSF) model, where effects of both 
finite interphase transport rates and turbulence 
on particle motion were considered using random 
sampling techniques. In general, the LHF and DSF 
models over- and under-estimated rates of particle 


spread and flow development, respectively. In 

■ contrast, the SSF model- yielded encouraging 
predictions of flow structure— except at high 
particle mass loadings, where effects of 
particles- on turbulence proper- ties (termed 
turbulence modulation by AlTaweel and 

-Landau?) which were not considered in the theory, 
were felt to be responsible for the deficiency. 

Relatively few theoretical studies treating 
particle-laden swirling flows have been published. 

■ Hamed 8 modeled the trajectories of solid particles 
in a flow downstream of swirler vanes. Dring and 
Suo 9 calculated particle trajectories in a free- 
vortex flow. Effects of turbulence on particle 
trajectories, however, were not considered in 
either of these studies. 

Due to the absence of existing measurements, 
the present experimental study was conducted to 
obtain data that could be used to assess multiphase 
flow models in a swirling flowfield. In order to 
obtain significant interaction between phases, 
relatively small particles (SMD = 39 ^m) were used. 
Emphasis was placed on obtaining gas-phase veloc- 
ities in the presence of particles without signal 
interference from particles. Initial conditions 
were measured near the jet exit. Single-phase 
jets were also studied to establish baseline 
results and are reported in Ref. 10. Also 
reported in Ref. 10 are predictions of particle- 
laden swirling jets based on a single particle 
size, the SMD. 

The present paper presents additional measure- 
ments and predictions for particle-laden swirling 
jets. Predictions are reported for the separated- 
flow models using seven particle size groups in 
order to better simulate the size distribution of 
the particles that were used in the study. Pre- 
dictions from the LHF model are not dependent on 
particle size; therefore, predictions from the 
LHF model reported in the present paper are iden- 
tical to those in Ref. 10. 


Experimental Methods 

Test Apparatus 

The particle-laden jets were directed ver- 
tically downward within a large enclosure (1.8 m 2 
by 2.4 m high) screened enclosure. The injector 
could be traversed vertically within the enclosure 
while the" enclosure and injector could be traversed 
together in the two horizontal directions. This 
arrangement allowed rigid mounting of all optical 
instrumentation used' during the study. 

The jet tube had an inside diameter of 19 mm 
and extended vertically downward for 100 injector 
diameters. Swirl was generated by injecting air 
tangentially through four 9.5 mm long slots, 
located 481.0 mm upstream of the jet exit. The 
swirl number (calculated as in Ref. 10) was changed 
by varying the amount of air injected through the 
tangential slots. Both the swirl and main air 
streams were metered with calibrated, critical- 
flow orifices. The solid-glass particles used 
during the study were injected into the flow far 
upstream of the injection tube using a vibrating, 
variable-speed screw feeder. The size distribu- 
tion of the particles used in the study is illus- 
trated in Fig. 1. Particle size generally ranged 
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from about 10 to 65 pm. Test conditions for the 
present study are summarized in Table I. Three 
particle-laden jets were studied with swirl num- 
bers ranging from 0 to 0.3. A single loading 
ratio of 0.2 was considered. 

Instrumentation 


Gas velocity . Gas-phase velocities were meas- 
ured in the presence of particles using a single- 
component phase/Ooppler particle analyzer, 
described in detail by Bachalo and Houser. 11 * 12 
With this instrument, particle size and velocity 
were measured simultaneously. Gas-phase velocity 
was measured by seeding the particle-laden jet and 
ambient surroundings with 1 ym nominal aluminum 
oxide particles and using the velocity measured 
for this size to represent the continuous phase. 
The green (514.5 nm) line of an Ar-Ion laser was 
used with a 602.4 mm focal-length lens. To reduce 
the size of the probe volume, a beam expansion 
ratio of 3 was used. The receiving optics were 
mounted 30° off-axis in the forward scatter direc- 
tion, with light collected using a 495 mm focal- 
length lens. Since frequency shifting was not 
available for this system; measurements of radial 
and angular velocities were not performed, and 
only mean and fluctuating axial velocities are 
reported . 

Particle velocity . Particle velocity was 
measured using two instruments. Particle veloc- 
ities, number averaged over all particle sizes, 
were measured using a two-channel (frequency- 
shifting on both channels) Ar-Ion laser velocime- 
try (LV) system. The receiving optics were placed 
~30° off-axis in the forward scatter direction. 

A 602.4 mm focal-length lens was used to collect 
the scattered light. Seeding particles were not 
introduced into the flow during the particle 
velocity measurements and both laser power and 
detector gain were reduced to insure that only 
signals from test particles were received. At 
each axial position, two traverses, 90° degrees 
apart yielded all required velocity components. 

Particle axial velocities were also measured 
using the phase/Doppler particle analyzer that was 
used to measure gas-phase velocities. For the 
particle velocity measurements, laser power was 
reduced to minimize signals from seeding material 
and to maximize signals from the particles. Mean 
and fluctuating axial velocities for particles 
ranging in size from 3.4 to 94.3 ym were measured. 

Particle mass flux . Particle mass flux was 
measured using an isokinetic sampling probe. 
Samples were collected on filter paper for a timed 
interval and weighed. Probes having Inside diam- 
eters of 2 and 5 mm were used to Insure adequate 
resolution and reasonable sampling times in var- 
ious regimes of the flow. Measured particle mass 
fluxes integrated across the jet were within +10 
percent of the calibrated particle flow rate at 
all axial locations. 

Theoretical Methods 


less than -0.5 to prevent any zones of recircula- 
tion. The boundary-layer approximations are 
adopted, however, the radial pressure gradient, 
which is usually neglected in the boundary-layer 
analysis, is considered. A k-e turbulence model 
is used to provide closure. The injector Mach 
number is <0.3; therefore, the kinetic energy and 
viscous dissipation of the mean flow are neglected 
with 1 ittle error. 

Three methods of treating multiphase flow, 
typical of current practice, are considered. The 
methods are: (1) locally homogeneous flow (LHF), 

where interphase transport rates are assumed to be 
infinitely fast and the flow can be treated like 
a single-phase, variable-density fluid; (2) deter- 
ministic separated flow (DSF), where finite inter- 
phase transport rates are considered but the 
dispersed phase is assumed to interact only with 
mean properties of the continuous phase; and (3) 
stochastic separated flow (SSF), where interphase 
transport rates and effects of turbulent disper- 
sion of the dispersed phase are treated. All 
three methods will be only briefly discussed since 
they have been fully described elsewhere. 1_ 6 

Continuous Phase 


Mean quantities for the continuous phase are 
found by solution of governing equations for con- 
servation of mass and momentum in conjunction with 
second-order turbulence model equations for tur- 
bulent kinetic energy and its rate of dissipation. 
The volume fraction of the particle phase was 
neglected, since void fractions (volume fraction 
of flow without particles) for the present flows 
exceeded 99.8 percent. The governing equations 
for the continuous phase can be put into the fol- 
lowing general form: 


h ' * ? I? <r P**> • r a, 




(1) 


The parameters 4> , S^, and S p ^, appearing in 
Eq. (1), as well as empirical constants, are sum- 
marized in Table II. 


A consequence of the presence of angular 
velocity is that even though the standard boundary- 
layer assumptions have been made, the radial 
momentum equation still must be considered: 

_~2 

_ £iL (2 ) 

3r “ r 


and cross-stream pressure gradients are not negli- 
gible. Due to the decay of angular velocity, the 
axial pressure gradient is also included. 

The turbulent viscosity was calculated as 
usual: 



(3) 


General Description 

The analysis is limited to steady, axisym- 
metric, dilute, sol id-particle-laden, weakly 
swirling turbulent jets in an infinite stagnant 
medium. The swirl number is restricted to values 


The flow leaving the injector was similar to 
fully developed flow and had no potential core. 
The boundary conditions for Eq. (1) are: 
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4 >*rw:r = 0, = 0; r -» », 4> = 0 


. (4) 


(9) 


= 0, rw = 0; r -» ®, rw = 0 


Initial conditions were measured at x/d = 0.5. 
e 0 was calculated from the definition of a tur- 
bulent length scale, as follows: 


Ck 3/2 



L 


(5) 


dt ' u pi 



Vp/ (Ul ’ 


V 


u - 


♦ a i 


(10) 


where i = 1,2,3 and the velocities shown in these 
equations are instantaneous velocities for a par- 
ticular eddy and particle group. Since the par- 
ticles are sol id. spheres , the drag coefficient was 
calculated as follows: 1 


where L was chosen as a fraction of the initial 
jet width to provide good agreement with the axial 
profile for k. Since only u and u' could be 
measured for the particle-laden cases, initial 
values of the continuous-phase angular velocity 
were estimated by subtracting the measured particle 
angular momentum from the single-phase values. 
Single-phase initial values of k and e were 
used for the particle-laden jet calculations. 

The present stochastic separated flow 
approach involves modification and extension of 
methods reported by Gosman and Ioannides. 12 
This involves computing particle motion as 
particles interact with a succession of turbulent 
eddies. Properties are assumed to be uniform 
within each eddy and to randomly change from one 
eddy to the next. At the start of a particle-eddy 
interaction, the velocity of an eddy is found by 
making a random selection from the probability 
density function (PDF) of veloci ty— assuming an 
isotropic Gaussian PDF having standard deviations 
(2k/3) 1/2 and mean values u, v, and W. A 
particle Is assumed to interact with an eddy for a 
time which is the minimum of either the eddy life- 
time or the time required for a particle to cross 
an eddy. These times are estimated following past 
practice, 2-5 assuming that the characteristic 
size of an eddy is the dissipation length scale: 


L 


e 


C 3/4 k 3/2 

P 


It 


( 6 ) 




Re < 1000; 


C D = 0.44, 


Re > 1 000 


Dispersed Phase 


( 11 ) 


SSF model . The dispersed phase was treated by 
solving Lagrangian equations for the trajectories 
of a statistically significant sample of individual 
particles (n groups defined by initial position, 
velocity and sample) as they move away from the 
injector and encounter a random distribution of 
turbulent eddies. This approach provides a means 
of treating effects of turbulent fluctuations on 
particle drag and dispersion as well as effects of 
particles on turbulence properties. Predictions 
presented in the present paper for the SSF model 
utilized 7000 particle groups. 

Key elements of the SSF model are the methods 
used to specify eddy properties and the time of 
interaction of a particle with a particular eddy. 

Particle source terms . The interaction between 
particles and the continuous phase yields source 
terms in the governing equations for conservation 
of axial and angular momentum. The source terms 
are found by computing the net change in momentum 
as each particle group 1 passes through computa- 
tional cell j 


and that the eddy lifetime is: 

t = L / ( 2k/3> 1 72 
e e 


(7) 




Z 


i = l 


n i m p 




( 12 ) 


Therefore, particles are assumed to interact with 
an eddy as long as both the time and distance of 
interaction satisfy the following criteria: 


At < Ax < L 
p - e p - e 


( 8 ) 



L 


i = l 



(13) 


Assumptions for particle trajectory calcula- 
tions are typical of analysis of dilute particle- 
laden flows: 1 drag is treated empirically, 
assuming quasisteady flow for spherical particles 
with no influence of nearby particles; particle 
collisions are neglected; effects of virtual mass, 
Bassett forces and Magnus forces are neglected 
with little error, since pp/p > 2000; and static- 
pressure gradients are negligible. Local ambient 
properties are fixed by Instantaneous eddy pro- 
perties, as described earlier, which implicitly 
provides for effects of turbulent fluctuations on 
particle dispersion and drag. 

With these assumptions, the position and 
velocity for each particle group can be found by 
integrating: 


where nj is the number of particles per unit 
time in each group. 

DSF model . Effects of turbulent fluctuations 
on particle drag and dispersion are ignored for 
the DSF model. Particle trajectories are found by 
integrating Eqs. (9) and (10), but the local mean 
velocity of the continuous phase replaces the 
instantaneous eddy velocity. Each initial condi- 
tion yields a single deterministic trajectory; 
therefore, 1400 particle groups suffice to numer- 
ically close the solution. Effects of particle 
drag in the mean momentum equations are found 
from Eqs. (12) and (13), similar to the SSF 
calculations. 
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LHF model . This approximation implies that 
both phases have the same instantaneous velocity 
at each point in the flow; therefore, the flow 
corresponds to a variable-density, single-phase 
fluid whose density changes due to changes in 
particle concentration. Turbulent dispersion of 
particles is then equivalent to that of a gas and 
particle inertia fully influences turbulence pro- 
perties; i.e., the method implicitly accounts for 
effects of turbulence modulation to the extent 
that the no-slip assumption is correct. 

The treatment of the variable-density fluid 
is similar to past practice, however f is defined 
here as the mass fraction of particles in the 
fluid. Measured initial values of f were used 
for the LHF predictions. Through the assumption 
of no-slip, there is no need to compute particle 
trajectories and all particle source terms in the 
governing equations for the continuous phase are 
zero. 

Numerical Solution 

The calculations for the continuous phase 
were performed using a modified version of GENMIX. 1 ^ 
The computational grid for the cases without swirl 
was similar to past work: 2-5 - 10 33 cross-stream 
grid nodes and streamwise step size was limited to 
5 percent of the current flow width or an entrain- 
ment increase of 5 percent-whichever was smaller. 

For the cases with swirl, 33 cross-stream grid 
nodes were also used, but streamwise step size was 
reduced to the smaller of either 2 percent of the 
current flow width or an entrainment increase of 
2 percent. The angular and radial momentum equa- 
tions required modifications to the standard 
solution procedure and this procedure is described 
in Ref. 10. 

Results and Discussion 

Predictions and measurements of the particle- 
laden jets are discussed in the following. As 
discussed earlier, only mean and fluctuating axial 
velocities could be measured for the gas phase; 
therefore, it was necessary to estimate initial 
values of angular velocity and turbulent kinetic 
energy for the gas phase. Initial angular veloc- 
ity of the continuous phase was estimated by sub- 
tracting the measured initial particle-phase 
angular momentum from the values obtained for 
single-phase cases. For both swirl flows, particle- 
phase angular momentum was -10 percent of the 
values measured for the single-phase cases. For 
the predictions shown here, initial values of k 
were assumed to be the same as the single-phase 
flows. Measured values of u' were ~20 percent 
lower across the entire jet width for the particle- 
laden cases (at the initial condition of x/d = 0.5) 
than for the corresponding single-phase flows. 
Predictions showed that a reduction of k of 
20 percent caused negligible changes in flow pro- 
perties except very close to the injector. 

Initial values of t for the particle-laden jet 
predictions were also unchanged from the single- 
phase cases. 

For the separated-flow models, the predictions 
utilized seven particle size groups ranging in 
size from 12.5 to 72.5 pm. Initial axial mean and 
fluctuating velocities for each particle size 
group were obtained from phase/Doppler measurements 
at x/d = 0.5. Initial mean and fluctuating 


angular and radial velocities were obtained from 
LV measurements at x/d = 0.5; thus each size 
group was given the same initial angular and 
radial velocity. 

For the LHF model, predictions are not depend- 
ent on particle size; therefore, predictions pre- 
sented in the present paper are identical to those 
reported in Ref. 10. 

Nonswirlinq Jet 

The predicted and measured variation of axial 
velocities in the streamwise direction for the 
nonswirling, particle-laden jet are illustrated in 
Fig. 2. Measured particle velocities reported in 
Fig. 2(b) are number averaged over all particle 
sizes. For comparison with these measurements, 
predictions from the separated-flow models were 
number averaged over all seven size groups as 
well. Because the particle loading ratio is 
relatively low, predictions from both the LHF and 
SSF models are nearly identical and show reasonably 
good agreement with the experimental measurements. 
For the gas-phase axial velocity, deterministic 
separated flow model predictions were identical 
to those from the stochastic separated flow model 
and are not shown. Particle axial velocity pre- 
dictions (Fig. 2(b)) using the LHF model and 
separated-flow models show distinct differences. 

The neglect of particle inertia, illustrated by 
predictions from the LHF model, overestimates the 
rate of decay of particle axial velocity. Predic- 
tions from both the stochastic and deterministic 
models are similar and show reasonably good agree- 
ment with measurements for particle axial velocity 
decay. 

Radial profiles of flow properties, number 
averaged over all particle sizes, for the particle- 
laden jets without swirl are illustrated in Figs. 3 
and 4. Measurements and predictions are 'illus- 
trated for x/d = 5 and 30. Gas-phase jet widths 
are slightly overestimated for the particle-laden 
flows at x/d = 30 (Fig. 4). The overestimation of 
gas-phase flow width is probably caused by turbu- 
lence modulation due to the particles. Even 
though the loading ratio was relatively low, values 
of u 1 for the particle-laden jets are reduced 
compared to the single-phase jets. 10 Particle 
axial velocities are reasonably well predicted by 
all the models, however, since the LHF model over- 
estimates the rate of axial velocity decay with 
streamwise distance, the values of axial velocity 
predicted by this approach are generally lower 
than the measurements. 

Radial profiles of particle mass flux 
(Figs. 3(c) and 4(c)) illustrate the different 
physical assumptions embodied in the three models. 

At x/d = 5, both separated flow models give similar 
predictions of particle mass flux. At larger 
streamwise distances, however, the neglect of tur- 
bulent dispersion of the particles causes the DSF 
model to underestimate the spread of the particles. 
Since mean gas-phase radial velocities are quite 
small at these axial distances, turbulent disper- 
sion is the only mechanism available for radial 
spread of the particles. The no-slip assumption 
of the LHF model causes the radial dispersion of 
the particles to be overestimated at x/d = 30. 

Only the SSF model, which accounts for both 
particle inertia and turbulent dispersion of the 
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particles, satisfactorily predicts the radial 
particle mass flux at x/d = 30. 

The predictions of fluctuating particle 
properties for the SSF model, also illustrated in 
Figs. 3(a) and 4(a), are reasonably good, however, 
u'p is underestimated. This is probably caused 
by the assumption of isotropic velocity fluctu- 
ations, which causes streamwise continuous-phase 
velocity changes experienced by the particle to be 
underestimated. In particular, measured values of 
u' are always greater than v' and w' for the 
single-phase jet 10 and it is expected that this 
behavior should be similar for the particle-laden 
jet. For the nonswirling jet, predictions using 
seven size groups are almost identical to those 
previously reported 10 using the SMD of the 
particle distribution. 

These results for the particle-laden jets 
without swirl are similar to those previously 
reported by other investigators for similar 
flows . 1- 6 

Measurements and predictions of particle 
mean axial velocity for 23, 43, and 63 pm par- 
ticles at x/d = 5, 15, and 30 are illustrated in 
Fig. 5 for the nonswirling jet. The data pre- 
sented were obtained using the phase/Doppler 
particle analyzer and represent axial velocities 
of particles from the size distribution shown in 
Fig. 1. Only predictions from SSF model are 
presented. As Illustrated in Fig. 5, decay of 
axial velocity of the larger particles is slower 
than the smaller particles. By x/d = 15, the 
63 pm particles have the highest velocity. Pre- 
dictions from the SSF model are In reasonably 
good agreement with measurements and correctly 
show the general trends. 

Swirl inq Jets 

Measurements and predictions for the 
particle-laden swirling jets are discussed In 
this section. Jets having swirl numbers of 
S = 0.16 and 0.3 were studied. 

Predicted and measured flow properties In the 
streamwise direction are illustrated in Figs. 6 to 
8 for the swirling, particle-laden jets. Pre- 
dicted and measured gas-phase axial velocities are 
shown In Fig. 6 for swirl numbers of 0.16 and 
0.3. As illustrated in Fig. 6, increasing swirl 
increases the rate of decay of axial velocity. 
Predictions from the LHF and SSF models are 
nearly identical and show reasonably good agree- 
ment with measurements for both swirling flows. 
Again, predictions of gas-phase axial velocity 
using the DSF model are nearly identical to 
predictions using the SSF model and are not shown 
in the figure. 

Particle mean axial and maximum angular 
velocities, number averaged over all particle 
sizes, are plotted as a function of streamwise 
distance in Figs. 7 and 8, respectively. As 
expected, neglecting slip between the phases 
causes predictions from the LHF model to over- 
estimate the rate of decay of particle velocities. 
Differences between predictions from the DSF and 
SSF models were small: both show reasonably good 

agreement with measurements. At the initial con- 
dition (x/d = 0.5), axial velocities of the 
particles were lower than the continuous phase 


except near the edge of the jet for both swirling 
flows. Particle axial velocities initially 
increase, due to momentum exchange from the con- 
tinuous phase, before beginning to decay. As 
i 1 1 ustrated in Fig. 7, both separated-flow models 
correctly predict this behavior. 

Radial profiles of flow properties for the 
particle-laden swirling jets are illustrated in 
Figs. 9 to 14. Velocity measurements, number 
averaged over all particle sizes, are illustrated 
in Figs. 9 to 12. Particle mean axial velocities 
for three size groups are illustrated in Figs. 13 
and 14 for the two swirling flows. 

Radial profiles of number-averaged flow pro- 
perties for the swirling, particle-laden jets at 
x/d = 5 are illustrated in Figs. 9 and 10 for 
S = 0.16 and 0.3 .respectively. As expected, 
increasing the swirl number increases the width 
of the particle-laden jet. The gas-phase flow 
width is slightly overestimated at x/d = 5 for 
both swirling flows. Again, this is probably due 
to turbulence modulation by the particles, as 
discussed earlier. Also, because the gas-phase 
flow width is overestimated, particle axial 
velocities are also overestimated for the 
separated-flow models. The SSF model over- 
estimates particle axial velocity to a greater 
extent than the DSF model. The LHF model, however, 
underestimates particle axial velocity, due to 
the neglect of particle inertia. At x/d = 5, 
particle axial velocity, number averaged over all 
size groups, is nearly equal to the continuous 
phase. For both swirling flows, the LHF model 
predicts that the maximum particle mass flux has 
shifted to the center of the jet, which is clearly 
not correct. The neglect of turbulent dispersion 
of the particles, illustrated by the DSF model, 
causes the particles to be confined to a rela- 
tively narrow region of the flow. If only mean 
properties of the continuous phase are con- 
sidered, the particles are transported by centrif- 
ugal forces to regions where radial velocity is 
small and tend to remain there. Considering tur- 
bulent fluctuations of the continuous phase gives 
better predictions of particle mass flus at 
x/d = 5. The predicted maximum particle mass flux 
is shifted radially outward when compared to the 
measurements. This is again probably caused by 
the overestimation of the jet width at this 
streamwise location. Predictions of mean particle 
angular and fluctuating velocities are quite good 
for the S = 0.16 flow, see Fig. 9(e). For the 
higher swirl number flow, S = 0.3, particle fluc- 
tuating axial velocities are overestimated (see 
Fig. 10), however, predictions of the other two 
fluctuating particle velocities show better agree- 
ment with measurements. Predicted particle fluc- 
tuating velocities are increased at all locations 
using seven particle sizes when compared to pre- 
dictions using the particle SMD, 10 especially 
for the higher swirl number flow. Predicted mean 
properties using the particle SMD are similar to 
those obtained during the present study using 
seven particle sizes. 

Radial profiles of gas-phase properties at 
x/d = 20 are illustrated in Figs. 11 and 12 for 
both swirl flows. Predictions are in better 
agreement with measurements at this position than 
closer to the injector, although the jet width is 
still si ightly overestimated. Predictions at 
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x/d = 20 are not as sensitive to initial condi- 
tions and the swirl component has almost com- 
pletely decayed. Again, there is little difference 
between predictions of gas-phase properties for 
the no-slip and separated-flow models. Predic- 
tions of axial velocity are in good agreement 
with measurements for all three models. However, 
since the LHF model overestimates the rate of 
decay of axial velocity, predicted unnormalized 
velocities from this model are lower than the 
measurements. Predicted values of u' p under- 
estimate the measurements, while v' p and w' p 
are in reasonably good agreement with measurements 
using the SSF model. Since effects of swirl have 
decayed at this axial location, ignoring the ani- 
sotropy of the continuous phase is the main reason 
for this behavior. Even at x/d = 20, predicted 
fluctuating particle velocities using a single 
particle size, the SMDjO are slightly lower 
than those of the present study using seven par- 
ticle sizes. 

The particle mass flux predictions 
(Figs. 11(c) and 12(c)) again highlight the dif- 
ferent physical assumptions of the three models. 
Particle mass flux measurements indicated that 
between x/d = 5 and x/d = 10, the maximum mass 
flux shifted to the center of the jet for both 
swirl cases. Since angular and radial velocities 
have decayed to relatively small values at this 
distance, and would tend to move particles out- 
ward, the only mechanisms for transport inward 
are turbulent dispersion and entrainment of 
ambient air at the edge of the jet. As shown in 
Figs. 11(c) and 12(c), the predicted maximum par- 
ticle mass flux from the SSF model has not com- 
pletely shifted to the center of the jet, however 
it is clearly evolving in this direction. In 
contrast, the DSF model predicts a very narrow 
distribution with no particles at the center of 
the jet. Compared to the nonswirling case, the 
LHF model underestimates particle dispersion for 
both swirling flows at x/d = 20. This behavior 
is caused by neglecting the radial and angular 
inertia of the particles which tend to transport 
them radial ly. 

Measurements and SSF model predictions of 
mean axial velocities for particle diameters of 
23, 43, and 63 pm in the two swirling flows are 
illustrated in Figs. 13 and 14. Again, measure- 
ments reported for each size group were obtained 
with the phase/Doppler particle analyzer. 

Radial profiles of u are shown at axial 
locations of x/d = 2, 5, 10, and 20. At 
x/d = 2, the velocity of the smaller particles is 
larger than the larger particles near the center 
of the jet. Because of their increased inertia, 
the axial velocity of the larger particles decays 
at a slower rate than smaller particles, and by 
x/d = 20, larger particles are moving at a higher 
velocity than the smaller particles at all radial 
locations. Increasing swirl Increases the varia- 
tion of velocity with particle size, see Figs. 13 
and 14. 

SSF model predictions of mean particle axial 
velocity for each size group show better agreement 
with measurements as distance from the injector 
exit increases. At streamwise distances of x/d 
= 10 and less, radial profiles of axial velocity 
were overestimated for both swirling flows. As 


discussed earlier, this is probably due to turbu- 
lence modulation since the continuous-phase jet 
width was also overestimated for these flows. 

At x/d = 5, particle velocities are under- 
estimated for both swirling flows (Figs. 13(b) 
and 14(b)). Since predicted centerline particle 
velocities were higher for the DSF model at 
x/d = 5 and 10, this appears to be the result of 
eddy specification in the stochastic model. This 
is especially true for the 23 pm particles. 

Eddies are assumed to travel at the mean gas-phase 
velocity at their initial radial position in the 
flowfield. Smaller particles, which have shorter 
relaxation times, tend to remain in a particular 
eddy longer than larger particles. 

Predictions from the SSF model display the 
correct trends with respect to variations in par- 
ticle velocity with size and streamwise distance. 
Predictions at x/d = 20 exhibit good agreement 
with measurements for both swirling flows. 

Conclusions 

Major conclusions concerning the models that 
were evaluated during this investigation are as 
follows: 

1. Predictions using the SSF model showed 
reasonable agreement with measurements for the 
particle-laden jets. In general, the LHF model 
overestimated the rate of particle velocity decay 
in the particle-laden jets due to the neglect of 
particle inertia. For nonswirling jets, the LHF 
model overestimated particle spreading rates due 
to the no-slip assumption. For swirling jets, 
the LHF model underestimated particle spreading 
rates due to the neglect of particle inertia. 

The DSF model underestimated particle spreading 
rates due to the neglect of the effect of turbu- 
lent fluctuations on particle motion. Only the 
SSF model, which accounts for both particle iner- 
tia and effects of turbulent fluctuations, cor- 
rectly predicted particle spreading rates over 
the entire flowfield. 

2. Particle axial fluctuating velocities 
were generally underestimated at far downstream 
axial locations. This is probably due to the 
assumption of isotropic velocity fluctuations in 
the SSF model, since fluctuating axial velocities 
are expected to be greater than fluctuating radial 
or angular velocities. 

3. Near the injector exit, jet widths were 
overestimated with the separated-flow models. 

This was probably caused by turbulence modulation 
by the particles, which was not considered in the 
analysis . 

4. Mean axial velocities for each particle 
size group were reasonably well predicted for the 
particle-laden jets using the SSF model. Predic- 
tions showed the same trends as measurements for 
particle streamwise axial velocity decay of each 
size group. 
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TABLE I. - SUMMARY OF TEST CONDITIONS 3 


Parameter 

Case 1 

Case 2 

Case 3 

Centerline air axial velocity, m/s* 3 

13.75 

11.9 

12.2 

Maximum air angular velocity, m/s c 

. 0 

2.84 

5.33 

Swirl number 

0 

0.16 

0.3 

Centerline particle axial velocity, m/s* 3 

13.8 

10.39 

10.26 

Maximum particle angular velocity, m/s* 3 

0 

1.48 

2.26 

SMD, pm d 

39 

39 

39 

Loading ratio e 

0.2 

0.2 

0.2 


aAmbient temperature and pressure, 296 K, 97 kPa; injector inside 
diameter, 19 mm; particle density, 2500 kg /m^. 

^Measured at x/d = 0.5. 

Estimated from single-phase cases. 

d The size distribution has a standard deviation of 15 pm. 
e Ratio of injected particle mass flow rate to air mass flow rate. 
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NUMBER 


TABLE II. - SOURCE TERMS IN GOVERNING EQUATIONS 




PARTICLE-LADEN JET 
S = 0 


U co = 13.75 M/S 
IT = 18 RT m/<j 



DIAMETER. (JM (b) PARTIC LE VELOCITY. 

FIGURE 1. - PARTICLE SIZE DISTRIBUTION. FIGURE 2. - STREAMWISE VARIATION OF NUMBER- 

AVERAGED AXIAL VELOCITY FOR THE PARTICLE- 
LADEN JET CS = 0). 
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u/u, 




r/x 

(d) PARTICLE FLUCTUATING VELOCITIES. 

FIGURE 3. - RADIAL VARIATION OF NUMBER-AVERAGED FLOW PROPERTIES AT X/d = 5 FOR THE 
PARTICLE-LADEN JET (S = 0). 
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(d) PARTICLE FLUCTUATING VELOCITIES. 

FIGURE A. - RADIAL VARIATION OF NUMBER-AVERAGED FLOW PROPERTIES AT x/d = 30 FOR THE PARTICLE-LADEN JET <S = 0) . 
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(b) S = 0.3. 

FIGURE 6. - STREAMWISE VARIATION OF GAS-PHASE 
AXIAL VELOCITY FOR THE PARTICLE-LADEN SWIRL- 
ING JETS. 



(b) S = 0.3. 

FIGURE 7. - STREAMWISE VARIATION OF NUMBER- 
AVERAGED PARTICLE AXIAL VELOCITY FOR THE 
PARTICLE-LADEN SWIRLING JETS. 
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(b) S = 0.3. 

FIGURE 8. - STREAMWISE VARIATION OF NUMBER- 
AVERAGED PARTICLE MAXIMUM ANGULAR VELOCITY 
FOR THE PARTICLE-LADEN SWIRLING JETS. 
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0 .1 .2 .3 .1) 

r/x 


(e) PARTICLE FLUCTUATING VELOCITIES. 

FIGURE 10. - RADIAL VARIATION OF NUMBER-AVERAGED FLOW PROPERTIES FOR THE 
SWIRLING (S = 0.3) PARTICLE-LADEN JET AT x/d - 5. 
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PARTICLE-LADEN JET 
S = 0.16 
x/d = 20 
u c = 3.05 m/s 
u pc = 3.693 M/s 

O DATA 

SSF 

DSF 

LHF 





(d) PARTICLE FLUCTUATING VELOCITIES. 

FIGURE 11. - RADIAL VARIATION OF NUMBER-AVERAGED FLOW PROPERTIES FOR THE SWIRLING 
(S = 0.16) PARTICLE-LADEN JET AT x/d = 20. 
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